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Abstract. This is a review on graphene quantum dots and their use as a host for spin 
qubits. We discuss the advantages but also the challenges to use graphene quantum 
dots for spin qubits as compared to the more standard materials like GaAs. We start 
with an overview of this young and fascinating field and will then discuss gate-tunable 
quantum dots in detail. We calculate the bound states for three different quantum 
dot architectures where a bulk gap allows for confinement via electrostatic fields: (i) 
graphene nanoribbons with armchair boundary, (ii) a disc in single-layer graphene, and 
(iii) a disc in bilayer graphene. In order for graphene quantum dots to be useful in the 
context of spin qubits, one needs to find reliable ways to break the valley-degeneracy. 
This is achieved here, either by a specific termination of graphene in (i) or in (ii) 
and (iii) by a magnetic field, without the need of a specific boundary. We further 
discuss how to manipulate spin in these quantum dots and explain the mechanism of 
spin decoherence and relaxation caused by spin-orbit interaction in combination with 
electron-phonon coupling, and by hyperfine interaction with the nuclear spin system. 



Quantum dots and spin qubits in graphene 



2 



1. Introduction 

Quantum dots are tiny islands of electrons where charging effects as well as quantum 
confinement play a crucial role [1]. Besides fundamental insights, these artificial atoms 
could also work as building blocks for the control of electronics at the single electron 
level. So called single-electron transistors are able to switch on and off electron transport 
through a dot by means of electrical gates using the effect of Coulomb blockade. The 
electron not only has a charge but also spin — its intrinsic angular momentum s = 1/2. 
The ability to use the spin rather than the charge of the electron to control electrical 
conduction is used in spintronics [2], a quite recent and very active field in mesoscopic 
physics. Exploiting the spin degree of freedom, a quantum dot can act as a spin filter 
[3, 4, 5] or as a spin-blockade device [6]. Spin being a quantum mechanical variable, can 
also be controlled in a coherent manner. Indeed, quantum dots have been proposed as 
a host for a quantum bit (qubit) where the spin of an electron residing in the quantum 
dot forms the fundamental two- level system of a qubit [7]. Arrays of such quantum 
dots with tunable tunnel-couplings between them would work as a universal quantum 
computer [8]. Recent progress in this field using GaAs-based two-dimensional electron 
gases is impressive [9] but some rather limiting sources of spin dccoherence have been 
identified and attributed to the rather strong spin-orbit interaction and to the hyperfine 
interaction with the nuclear-spin system of the host material. 

In this review, we would like to give an overview on the field of quantum dots in 
graphene — a single sheet of graphite [10, 11]. The first part of the review focuses on our 
recent theoretical work on quantum dots in single-layer and bilayer graphene that can be 
controlled and formed by electrical gates. These proposals solve two main problems of 
graphene quantum dots as building blocks for electron spin-based quantum computing: 
The challenge of achieving a controllable confinement and the reliable breaking of the 
valley degeneracy. The second part of the review is devoted to spin manipulation 
and decohcrence in graphene quantum dots. Here, we put an emphasis on the 
two most prominent spin decoherence channels: spin-orbit interaction in combination 
with electron-phonon coupling to lattice vibrations and hyperfine interaction with the 
surrounding nuclei. 

The review is organized as follows: We start with a general overview of the field 
of graphene quantum dots in Section 2. In Section 3, we will discuss the bound states 
as a function of confinement and/or as a function of magnetic field in three different 
geometries: (i) graphene nanoribbons, (ii) discs in single-layer graphene, and (iii) discs in 
bilayer graphene. The three proposals (i)-(iii) have in common that electrical gates can 
define the spatial extend of the quantum dot, however the way the valley degeneracy of 
graphene is broken — a prerequisite for spin qubits in graphene — is different for the three 
proposals. In Section 4, we present and discuss the physics of spin qubits in graphene 
quantum dots including an extended discussion of spin decoherence properties. Finally, 
in Section 5, we conclude and give a brief outlook. 
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2. General overview of graphene quantum dots 

Graphene is the first truly two-dimensional condensed matter system discovered 
experimentally in the group of A. Geim at the University of Manchester in 2004 [12]. 
The low-energy dynamics of electrons is governed by states near two inequivalent K- 
points in the Brillouin zone — so called valleys — and resembles the one of massless Dirac 
fermions [13, 14, 15, 16] where a new pseudospin degree of freedom appears due to the 
two sublattices defining the honeycomb lattice of graphene. This low-energy Dirac 
Hamiltonian indeed describes the physics of electrons and holes in the vicinity of the 
Dirac points as experimentally confirmed by the existence of the half-integer quantum 
Hall effect in graphene [17, 18]. 

Quantum dots in graphene have been extensively discussed in the recent literature 
- both from the theory side as well as from the experimental side. Theoretically, many 
approaches to graphene quantum dots exist: Graphene islands can be created by cutting 
flakes of graphene into the desired shape. In this case, boundary conditions and chaotic 
behavior become relevant [19, 20]. On the other hand, alternative approaches exist 
for electrostatically defined quantum dots, including quasi bound states in single-layer 
[21, 22, 23, 24] and bilayer graphene [23] as well as true bound states in single-layer 
[25, 26, 27] and bilayer graphene [28, 26, 29]. An electrostatic potential alone cannot 
produce true bound states unless some very specific symmetries arc met [24] because 
of the phenomenon of Klein tunneling [30, 31]. To overcome this problem, one needs 
to find a way to generate a gap in the bulk of graphene. Another possibility to confine 
Dirac fermions is by using spatially inhomogeneous magnetic fields [32]. Furthermore, 
quantum antidots, that can be created in graphene by cutting out holes, have been 
proposed as potential host systems for spin qubits [33]. 

Experimentally, graphene quantum dots are created by etching or scratching 
graphene islands attached to leads via constrictions. Single electron transistor behavior 
[34, 35] as well as quantum confinement [35, 36] have been observed, more recently also 
in double quantum dots [37, 38]. Energy levels have also been studied in magnetic fields 
where indication of the formation of Landau levels was observed [36, 39]. The zeroth 
Landau level then can be used to locate the electron-hole crossover in graphene quantum 
dots [39]. Observation of quantum confinement as well as Coulomb blockade effects in 
electron transport requires that an electron is trapped in the dot long enough such 
that size quantization and charging energy effects become relevant. The constrictions 
separating the dot from the open contacts act as tunnel barriers. It is believed that the 
constrictions are tiny nanoribbons where a gap is formed in graphene due to confinement. 
Graphene nanoribbons in combination with a back gate and top gates themselves can 
be used to form single quantum dots [40] as well as double quantum dots [41]. To 
understand the behavior of the level structure in edged graphene quantum dots and 
graphene nanoribbons, edge disorder as well as Coulomb blockade effects play a crucial 
role [42, 43, 44]. Due to the large disorder induced by the edges it is rather difficult to 
control the level-structure of the quantum dots. Especially in graphene nanoribbons. 
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disorder plays a dominant role [45, 46]. Nevertheless, new experiments report the 
observation of excited states [36] as well as spin states [47], pointing towards some 
improvements in control and stability of graphene quantum dots. 

We mention in passing that it has been theoretically shown that edge states, 
prominent in graphene nanoribbons with zigzag edges, are present in any quantum 
dot shape with sharp boundaries [48, 49] which seemingly predicts a rather universal 
feature of graphene nanoflakes. 

3. Electrostatically defined quantum dots in graphene 

In this section, we discuss our approach to use electric fields to confine electrons in 
graphene and bilayer graphene. A straightforward approach as used in GaAs quantum 
dots [9] is not possible due to Klein tunneling. Several ways are possible to induce a 
gap in bulk graphene. In general, quantum confinement can lead to the opening of 
a gap in ribbons [50, 51, 52]. Within the tight- binding approximation of graphene, 
armchair boundary conditions can lead to an insulator and gate-tunable quantum dots 
in ribbons can be created. Another promising direction is to start with bulk graphene 
and induce a gap via the interaction with a substrate [53, 54, 55]. In the following, we 
discuss three quantum dot architectures that allow for bound states in graphene tunable 
by electrostatic fields: (i) graphene nanoribbons with armchair-terminated boundaries, 
(ii) discs in single-layer graphene, and (iii) discs in bilayer graphene. Special emphasis 
is given on the ability to controUably break the valley degeneracy, a prerequisite for 
two-qubit gates [25, 26] for spin-based quantum computing in graphene. 

3.1. Quantum dots in graphene nanoribbons 

We first concentrate on a single quantum dot which is assumed to be rectangular with 
width W and length L, see Fig. 1. The basic idea of forming the dot is to take a ribbon 
of graphene with semiconducting armchair boundary conditions in x-direction and to 
electrically confine particles in y-direction. 

The low energy properties of electrons (with energy e with respect to the Dirac 
point) in such a setup are described by the 4x4 Dirac equation 

I y -cr^Ox + Oydy J 

where the electric gate potential is assumed to vary stepwise, = //gate in the dot 
region (where < y < L), and //(y) = /^barrier in the barrier region (where y < or 
y > L). In Eq. (1), and ay are Pauli matrices (denoting the sublattices in graphene). 
The four component spinor envelope wave function ^ = (^^-*, — \ — -*) 
varies on scales large compared to the lattice spacing. Here, A and B refer to the 
two sublattices in the two-dimensional honeycomb lattice of carbon atoms, whereas 
K and K' refer to the vectors in reciprocal space corresponding to the two valleys in 
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Figure 1. Quantum dot in graphene nanoribbon. A ribbon of graphene with 
semi-conducting armchair boundaries is schematically shown. Two barrier gates (blue) 
define the rectangular size of the quantum dot (with width W and length L). A back 
gate (red) allows one to shift the energy levels in the dot. Two or more quantum dots 
of this type can be easily put in series in a single nanoribbon as will be discussed in 
more detail below (after Ref. [56]). 



the bandstructure of graphene. The appropriate semiconducting armchair boundary 
conditions for such a wave function can be written as [a = A, B) [51] 

^L''\=w = e^'-/'¥Pl=^. (2) 

These boundary conditions couple the two valleys and are, thus, the reason why the 
valley degeneracy is lifted [57]. It is well known that the boundary condition (2) yields 
the following quantization conditions for the wave vector = Qn in x-direction [51, 58] 

g„ = (n ± l/3)7r/W, n e Z. (3) 

The level spacing of the modes (3) can be estimated as Ae ~ hvir/SW, which gives 
Ae ~ 30meV, where we used that v ~ lO^m/s and assumed a quantum dot width of 
about W ^ 30 nm. Note that Eq. (3) also determines the energy gap for excitations as 
Egap = ^hvqo- Therefore, this gap is of the order of 60 meV, which is unusually small for 
semiconductors. This is a specific feature of graphene that will allow for long-distance 
coupling of spin qubits as will be discussed in Sec. 4 below. 

We now present in more detail the ground-state solutions, i.e. n = in Eq. (3). The 
corresponding ground-state energy e can be expressed relative to the potential barrier 
= yUbarricr in the regions y < and y > L as e = /^barrier ± ^v{qQ + fc^)^/^. Here, the ± 
sign refers to a conduction band (+) and a valence band (— ) solution to Eq. (1). For 
bound states to exist and to decay at y ±oo, we require that hvqo > \e — /ibarrierl, 
which implies that the wave vector ky = k in ^/-direction, given by 



k 



- ((^ - yUbarricr)/^^^)^ (4) 
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Figure 2. Quantum dot levels. Bound-state solutions of a dot witli aspect ratio 
qqL = ttL/3W = 5. The diagonal lines indicate the region in which bound-state 
solutions do exist given by Eq. (5). All energies are taken in units of hvqo (after 
Ref. [56]). 



is purely imaginary. In the dot region {0 < y < L), the wave vector k in ?/-direction is 
replaced by k, satisfying e = /igate±^'i'(5'o + ^^)^''^- Again the ± sign refers to conduction 
and valence band solutions. (In the following, we focus on conduction band solutions 
to the problem.) In the energy window 

arrier ; 

(5) 

the bound state energies are given by the solutions of the transcendental equation 



, frT\- h'vk^/ihvqoy -{e- /ibarrier)^ 

ian[KL — — —. [0) 

[£ - /Ubarricr)(£^ - /igate) " (/t^go) 

We show a set of solutions to Eq. (6) for a dot with aspect ratio q^L = nL/SW = 5 in 
Fig. 2. 



3.2. Graphene disc in single-layer graphene 

In this subsection, we discuss a quantum dot in gapped graphene defined via electrostatic 
gates without requiring a certain boundary condition as the quantum dot will be purely 
defined electrostatically by a potential U{x,y). The interaction with a substrate can 
induce different potentials for the two sublattices thereby breaking inversion symmetry 
[53, 54, 55]. This interaction gives rise to a mass term Maz in the Dirac equation. 
Then the Hamiltonian in the valley-isotropic form [59, 60] and with the inclusion of a 
perpendicular homogeneous magnetic field is [26] 

H^ = Ho + TAaz + U{x,y), (7) 

where = v{p + eA) ■ cr, B = V x A = (0,0,5), v = 10^ m/s is the Fermi 
velocity and r = ± differentiates the two valleys K and K'. We choose the symmetric 
gauge A = y(— ?/, x,0) and assume a circular symmetry in the confinement potential 
U{x,y) = U{r) with r = y^xM-~y^. In the following we set h = l. 
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Figure 3. Left: Energy diagram for quantum dot in single-layer graphene. A 

constant mass term A and an electrostatic potential with height Uq give rise to bound 
states (dashed line) in the conduction band (c) defining a quantum dot of radius R. 
Note that the confining potential U{r) is repulsive for holes in the valence band (v). 
Right: Bound state levels as function of dot radius. Energy levels for Uq — A 
and for j — 1/2 at zero magnetic field. Full lines correspond to r = +1, dashed lines 
to T = -1 (after Ref. [26]). 



Since Hr commutes with the total angular momentum operator = —id^, + <7z/2, 
the energy eigenspinors can be chosen to be eigenstates of 

with j the eigenvalue of Jz which has to be an half-odd integer. To find the bound 
states of Eq. (7) we solve the radial Dirac equation HT-{r)x^{r) = Ex'^{r), with 
X'^ir) = iXAir),XBir)V and 



+ br 



H,(r) = -iva.dr + rAa, + f/(r) + vay\ ^ ^ + for ' ' 

where h = eB/2. The radial-part of the spinor solutions for bound states inside 
{U{r) = 0) and outside {U{r) = Uq) the quantum dot are expressed in terms of 
hypergeometric functions M{a,b,z) (r < R) and U{a,b,z) (r > R) [26], where R is 
the dot radius, see Fig. 3. 

Matching the solutions at r = i?, we find the characteristic equation for the energy 
eigenvalues E of quantum dot bound states. For j > 0, we obtain 

e+M(g<,j + 1/2, x) f/(g>,j + 3/2, x) 

-e^M(g<,j + 3/2, x) [/(g>,j + 1/2, x) = 0, (10) 

and for j < we obtain 

e- M(g< , -J + 3/2, x) U{q^ - 1, -J + 1/2, x) 
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-e< M(g< - 1, -j + 1/2, x) U{q>, -j + 3/2, x) = 0, (11) 

where x = hB?—{l/2){R/lBY with = sJh/eB the magnetic length. Without loss 
of generality, we choose here B positive. The bound state levels for negative B are 
obtained from the symmetry Hr{i, B) = H^r{—j, —B). We have further introduced the 
parameters g<,> = (j - 1/2) ^(j) + 1 - (e|^> - A'')/Abv\ i+ = (e< - rA)/4(j + 1/2), 
e+ = 6t;V(e> + rA), e< = (j - l/2)/(e< + rA) and e> = l/{e> + tA) with e{x) the 
Heaviside step function, and = E, = E — Uq. 

In Fig. 3 (right panel), we show the evolution of quantum dot levels as a function of 
dot size for angular-momentum quantum number j = 1/2. Pull lines and dashed lines 
correspond to the two valleys. Due to the symmetry E{j,T) = E{—j, — r), the two set 
of curves display also the cases j — 1/2 and j — —1/2 in the same valley. The different 
solutions for the dashed and full hues are therefore a direct consequence of effective 
time-reversal symmetry (eTRS) breaking in a single valley at zero magnetic field by a 
finite mass term [61], sec Sec. 3.4. 

In Fig. 4, we show the bound states of the quantum dot as a function of magnetic 
field evaluating the characteristic equations Eqs. (10) and (11) numerically. In Fig. 4(a), 
we show the low-lying bound states in the conduction band. Note that the valley- 
degeneracy (or orbital degeneracy) is broken at finite magnetic field. The largest level 
spacing between the (non-degenerate) groundstate and first excited state we estimate 
from Fig. 4(a) to be at R/Ib ~ 1-8 and is about 165 meV/it![nm] for the parameters 
used in Fig. 4. At R/Ib ~ 1-8, we obtain for the Zeeman splitting A^, = g/J-sB ~ 200 
meV/i?^[nm] using g — 2 which shows that the level spacing is always larger than the 
Zeeman energy for reasonable dot sizes. Considering a QD with R = 25 nm, we obtain 
a valley splitting A^k' at R/Ib ~ 1-8 of about 6.6 meV corresponding to 77 K, being 
much larger than 4 K, the temperature achieved by cooling with liquid helium. The 
necessary magnetic field corresponding to R/Ib = 1-8 is 5=3.41 T (and B = 0.85 T 
for i? = 50 nm with Akk' ~ 3.3 meV) which is also easily achievable in the laboratory. 
A gap of size 0.23 eV has been concluded from ARPES data in graphene on top of a 
SiC substrate [54]. Therefore, the gap A and also the confining potential step height 
Uq could easily be larger than the QD level spacing S = hv/R which is about 26 meV. 
These results suggest that such QDs confined in graphene could be an ideal host for spin 
qubits where the orbital degeneracy is controllable by a magnetic field. Note that these 
quantum dots will not suffer from edge roughness as much as quantum dots defined in 
graphene nanoribbons. 

In Pig. 4(b), we show the merging of the quantum dot states with the bulk Landau 
levels (LLs) 

En = ±5^{A/5f + 2n{R/lBf] n = 1, 2, 3, ... (12) 

with increasing magnetic field. Note in particular, that there is a zero mode LL at 
E = — rA which lies entirely in one valley [16]. In contrast, in Ref. [27], the zero mode 
LL is at = as a result of the assumption of a finite mass outside the dot, but zero 
mass inside the dot. 
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Figure 4. Dot levels at finite magnetic fields, a) Numerical evaluation of 
characteristic equations (10) and (11) as a function of R/Ib with Ib = {h/eBy^'^ 
the magnetic length and R the dot radius. We use A = 10(5 with S — hv/R and 
Uq = A. a) At small i?-fields we observe a breaking of the level degeneracy, shown for 
(j, r) = (i,±l), (— i,±l), and (±|,±1). The full lines are for r = 1 and dashed lines 
are for r = — 1 corresponding to the two valleys of graphene. b) Same parameters as in 
a), but for larger magnetic fields. The energy levels converge to the bulk Landau levels 
with increasing R/Ib- Included are levels for j ^ zt-. 



±|, and ±1 (after Ref. [26] 



3.3. Graphene disc in bilayer graphene 

Bilayer graphene is the two-layer analog of the single layer, which in addition is coupled 
by a tunneling matrix element t± coupling two sublattices of different layers (the so- 
called Bernal stacking). A voltage V between the two layers breaks inversion symmetry 
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(like the mass term A in the single layer) and opens a gap proportional to the voltage 
[62, 63]. In addition, the combination of a top gate and a back gate allows to tune the 
gap and the average potential U{r) independently, see Fig. 5. 



bilayer graphene 




dopant atoms 
.,/,.. 



SiOo 



back gate 



Figure 5. Quantum dot in bilayer graphene. A back gate and dopants on top 
of the bilayer control the voltage V between the layers — leading to a controllable gap 
opening — as well as the Fermi energy (band filling). An additional top gate allows to 
induce a spatially inhomogeneous electrostatic potential U{r) analogous to the single- 
layer model which leads to bound states in the conduction band (or valence band) of 
the bilayer. Another possibility is to use a split top gate (instead of a combination of 
top gate and dopants) to achieve a similar confinement (after Ref. [26]). 



To investigate gate-tunable quantum dots in bilayer graphene, we separate the 
Hamiltonian in the bilayer into two parts: H = Tio + Til- encodes the motion of 
the electrons within the planes and is given by two copies of the Dirac equation. In the 
valley-isotropic representation it takes on the form {h = v = 1) 



Ho 



( 



\ 





Vx - iPy 






+ iPy 
















Px - iPy 





(13) 



Px + iPy 

in both valleys. Like in the case of the single layer we add a magnetic field by the 
minimal coupling prescription p — (p + eA) with A = {B /2){—y, x, 0). The other part 
of the Hamiltonian (i.e. HI) encodes the biasing field and the hopping t± between the 
two planes. The interplane hopping matrix element tj_ has recently been measured to 
be t± = 0.40 eV [64, 65]. In the simplest approximation, we may take 

\ 








tV 
2 





t± 







tV 
2 







t± 





tV 
2 


V 












tV 
2 



+ f/(r)l. 



(14) 



with f/(r) the applied electrostatic potential profile as defined in the previous section. 
The index r = ±1 again distinguishes the two valleys (note that in the valley-isotropic 
representation the basis is chosen such that the two planes in the bilayer are exchanged in 
the spinors that describe different valleys) . The same Hamiltonian has been analyzed 
in Ref. [28] at zero magnetic field and in Ref. [29] at finite magnetic field. However, 
the confinement described in Eq. (14) by U{r) has been achieved in Refs. [28, 29] by a 
position dependent "mass term" V{r), instead, leading to qualitative and quantitative 
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differences as the gap is absent in the dot (i.e. V = in the dot center), and electrons 
as well as holes can be confined in the dot, simultaneously 

To diagonalize "H (i.e. to find the eigenspinors \E' that fulfill 'H\I' = E^^) we go 
to cylindrical coordinates in which the states are easily classified according to their 
conserved value of total angular momentum m (m being an integer). More explicitly, 
we factor out the angular dependence of the states according to 





f 


1 











\ 





























1 






















) 



(15) 



Note that the angular momentum in the bilayer case is an integer m, in contrast to 
the half-odd integer j in the single layer case, which reflects the different pseudospins 
in the bilayer (pseudospin 1) and single- layer (pseudospin 1/2). With the definitions 
j — m-\-l/2 and s — sgn(S), the Hamiltonian I-Lqi which now acts on ^i, can be written 
as 

1 



X 



/ 


V 









si 









dp 







i 
i 




si 



(16) 



/ 



In the latter equation, we have introduced the dimensionless coordinate ^ = r/(\/2lB), 
where Ib = \Jfil (e|-B|) is the magnetic length. The eigenvalue problem can now be 
solved by using the general solutions of the ordinary differential equation imposed by "Ho- 
The details of finding the eigenvectors and eigenenergies are somewhat more involved as 
in the single layer case, and we refer the interested reader to Ref . [26] . The eigenspinors 
are again expressed in terms of hypergeometric functions M(a, 6, z) (inside the dot) and 
C/(a, 6, z) (outside the dot). For a fixed energy two solutions for constant U{t) exist 
(corresponding to the two bands of the bilayer at positive £'), so that the bound-state 
solutions are found by matching linear combinations of pairs of solutions inside and 
outside the dot at r = i? which requires to find the determinant of a 4x4 matrix. 

The results for the bound state energies as a function of magnetic field are shown 
in Figs. 6 and 7. The most important result of our study can be seen in Fig. 6 where 
we display the energy levels of a dot as a function of the magnetic field. At zero 
magnetic field, the degeneracy of the levels in the two valleys is clearly displayed. With 
increasing the magnetic field, the orbital degeneracy is lifted. The symmetry of the 
levels is analogous to the case of the single layer discussed above. The states that are 
degenerate at zero field are related by time-reversal symmetry which means that they 
correspond to opposite values of angular momentum ±m in different valleys. The typical 
eTRS of ±m within one valley is already broken by a finite "mass term", similar to the 
case of a disc in single- layer graphene discussed in Subsection 3.2. 
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Figure 6. Energy levels of bilayer quantum dots. Energy levels in a relatively 
small bilayer quantum dot (radius i? = 25 nm) as a function of the magnetic field. 
Parameters are as follows: t± = 0.4 cV = 15.19hv/ R, Uq = 1.52hv/R and s = 1 
(i.e. positive B-field). The solid and dashed lines are for different valleys, a) For 
V = IMv/R, b) for V = 6hv/R (after Ref. [26]). 



An important feature of the bilayer as opposed to the single layer is the 
unconventional Mexican hat-like dispersion relation near the band edge. This is most 
apparent for a large value of the bias field V, the case shown in Fig. 6b). It is clear 
that there are many closely spaced levels near the band edge. This is a feature of the 
enhanced density of states near this particular energy [66]. 

For a large quantum dot, it is also possible to reach the regime where the dot levels 
are described by the Landau levels. This feature is seen in Fig. 7a) where we display 
the bound states for m = 0,m = ±1 for large magnetic fields. Note that the dot levels 
tend to approach the bulk Landau levels displayed in Fig. 7b). 

3.4- Symmetries and breaking of orbital degeneracy in graphene quantum dots 

In this subsection, we discuss the symmetries of graphene and their breaking by a 
smooth boundary which is the case for quantum dots created by electrical gates, see 
Subsections 3.2 and 3.3. Here, we specialize to the case of single-layer graphene (similar 
arguments hold for the bilayer). The main point of Subsections 3.2 and 3.3 has been to 
demonstrate that the valley- degeneracy can be lifted without mixing the valleys. In the 
following, we provide the reason using general symmetry arguments. The Hamiltonian 
Eq. (7) breaks two symmetries present in pure graphene (i) the effective time-reversal 
symmetry p — )■ — p and a — t- —a, (ii) inversion symmetry x — t- — x of the graphene 
lattice. Both symmetries (i) and (ii) would lead to degeneracies in the spectrum — 
(i) within the same valley, (ii) within different valleys. In the case, where the valleys 
remain uncoupled by the boundary of the dot, there is another fundamental symmetry 
that leads to a degeneracy of states in different valleys, which is time-reversal symmetry 
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Figure 7. Dot levels at high magnetic fields. Merging of bilayer quantum dot 
levels to the bulk LLs as function of magnetic field for a relatively large bilayer QD with 
R=67.48 nm and tj_ = 0.4eV = ilhv/R, Uq = 3.5hv/R, V = h.lihv/R for m = 0, ±1 
and s = 1 (i.e. positive i?- field), b) Bulk Landau Levels (LL) which are approached 
by the dot levels in a) almost perfectly at high fields in this parameter regime. Full 
lines are for r = +1 and dashed lines are for r = — 1 (after Ref. [26]). 



T = —{Ty® <yy)C [59]. However, time- reversal is easily broken by applying a magnetic 
field. Effective time-reversal symmetry (eTRS) T = icjyC and inversion symmetry 
X = —{t.j; ® az)!^ are both broken by the mass term Ar^ ® Oz-, where C and Ti, denote 
complex conjugation and space inversion, respectively. 

The magnetic field in z-direction alone breaks T and T, but not X. Only the 
combination of a mass-term (generated by the substrate in single-layer graphene and 
the voltage V between the two layers in bilayer graphene) and the magnetic field breaks 
all the orbital degeneracies. In Subsections 3.2 and 3.3, we have provided two realistic 
systems where these symmetries are broken. 

Additionally, we mention here that in single-layer graphene the symmetries T and X 
are generally broken by any termination of the graphene lattice that does not couple the 
two valleys. At the boundary, the spinor needs to satisfy the linear equation i\} = Aiip 
where Ai is in general a 4 x 4 matrix of the form Ai = z/ ■ r ® ■ a where z/ is a 
vector on the Bloch sphere and is a vector in the plane tangential to the boundary 
[48]. One easily shows that [M, T] ^ and [A^,X] 7^ if the vector u = ±z (i.e. if the 
valleys are not coupled). Then, it is sufficient to additionally break T by an external 
magnetic field if one wants to lift the valley-degeneracy [67]. 

4. Spin qubits in graphene 

Since the pioneering proposal on spin qubits for solid state quantum information 
processing [7], astonishing experimental achievements have been made, particularly 
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in GaAs-based quantum dots, see Ref. [9] for a recent review. In these (up to now) 
most advanced spin qubit systems, two major sources of spin decoherence have been 
identified: (i) spin-orbit interaction in combination with electron-phonon couphng to 
lattice vibrations, and (ii) hyperfine interaction of the electron spin with the surrounding 
nuclear spins. 

Generally speaking, the ideal qubit should be easy to manipulate and should couple 
rather little to its environment. Therefore, it is natural to think about new host materials 
where at least one of these two points can be optimized. We show below that spin 
qubits in graphene quantum dots gain on both sides. They have more flexibility on the 
manipulation side and are superior as far as spin decoherence is concerned as compared 
to more conventional semiconductors like GaAs. 

4.1. Manipulation of spin qubits in graphene quantum dots 

For universal quantum computing, singlc-qubit and two-qubit manipulations are 
necessary. Single qubit rotations of spin qubits are naturally done by electron spin 
resonance (ESR) - as proposed in Ref. [68] and measured in Ref. [69] - or by electric- 
dipole-induced spin resonance (EDSR) - as proposed in Ref. [70] and measured in 
Ref. [71]. The Rabi frequency /Rabi at which the qubit rotates, for instance, in the 
ESR experiment [69] is proportional to the electron spin gf-factor /Rabi = g/J'sBac/'^h 
where fj,B is the Bohr magneton and B^c the external oscillating magnetic field used to 
rotate the spin. Notably, the electron spin g^-factor differs for different materials. In 
GaAs quantum dots, it has been measured to be \g\ < 0.43 [72] whereas, in graphene 
quantum dots, it has recently been determined to be close to \g\ = 2 [47]. Thus, it 
is possible to rotate the electron spin in graphene quantum dots using ESR about five 
times faster than in GaAs quantum dots using the same field strength of the external 
oscillating magnetic field. This is an important gain because all qubit manipulations 
need to be done fast to be able to implement fault-tolerant quantum computing [73]. 

Another important advantage of graphene spin qubits is related to the small band 
gap in graphene nanoribbons. (For a ribbon width of about 30nm, the band gap can be 
estimated to be of the order of 60 meV.) This fact yields additional flexibihty for two- 
qubit operations. Two-qubit operations are usually done via the Heisenberg exchange 
interaction iJcx = -/Si -82 where Sj (z = 1, 2) denotes the spin of the coupled qubits [74]. 
Within a Hubbard model approximation with t ^ f/, where t is the tunneling matrix 
element between the quantum dots and U the onsite Coulomb energy, the singlet-triplet 
splitting J can be estimated to be J At^ /U . The value of U is usually fixed by the 
device geometry and the electromagnetic environment of the quantum dot. For graphene 
quantum dots, typical values of U are about 10 meV, see, for instance, Ref. [45]. The 
tunneling matrix element t can, however, be easily tuned by increasing or decreasing 
the overlap of the wave functions of the electrons in the two quantum dots. In graphene 
or better to say in small band gap semiconductors, this manipulation can be done in 
two distinct ways. Either through tunneling via conduction band states (i.e. normal 
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Figure 8. Long distance coupling of three graphene qubits. The energy bands 
of three quantum dots (including barrier regions) are shown in which dot 1 and dot 3 
are strongly coupled via cotunnelling processes through the valence bands of barrier 2, 
barrier 3, and dot 2. (The valence bands are shown in red and the conduction bands 
in blue.) Importantly, the center dot 2 is decoupled by detuning. The energy levels 
are chosen such that Ae2 <C Aei (after Ref [25]). 

tunneling) or through tunnehng via valence band states (i.e. Klein tunneling). This 
has been predicted for graphene nanoribbons in Ref. [25] and experimentally realized in 
carbon nanotube quantum dots in Ref. [75]. 

The most important physical consequence of this additional flexibility is the 
appearance of a new type of long-distance coupling between graphene spin qubits as 
illustrated in Fig. 8. By means of Klein tunneling, two distant qubits can be strongly 
coupled without touching the states of intermediate qubits that might be located 
between the two. Thus, a ribbon of graphene hosting many spin qubits in a line can be 
viewed as a qubit piano where any two of them can be coupled with leaving the states of 
the others unchanged, see Fig. 9 for a schematic illustration. Interestingly, this feature, 
i.e. the availability of non-local interactions, is important for quantum error correction 
since it raises the threshold for fault-tolerant quantum computing [77]. 

In the previous section, we have given convincing arguments why the spin qubit 
manipulation in graphene has decisive advantages as compared to other commonly used 
host materials. In the next section, this advantage will be complemented by promising 
spin relaxation and dephasing properties of electrons and holes in graphene quantum 
dots. 

4-2. Spin relaxation and dephasing in graphene quantum dots 

Why can we expect stable spin qubits in graphene quantum dots? There is hope that 
spin relaxation and dephasing will be very weak in graphene for the following reasons: (i) 
Carbon is a light element with atomic number 6. Hence, its atomic spin-orbit interaction 
is weak as compared to heavier elements. However, such a statement should be taken 
with care because, in the solid state, spin-orbit coupling is often times dominated by 
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Figure 9. Qubit piano. Illustration of many spin qubits in a line hosted within a 
graphene nanoribbon. Quantum dots arc red bars and barrier regions are blue bars. 
Different spin qubits that are strongly coupled to each other via Klein tunneling are 
marked with the same color (after Ref. [76]). 

bulk inversion or structure inversion asymmetry. Therefore, crystal structures of light 
elements can (under certain circumstances) exhibit rather strong spin-orbit coupling. 
Prime examples are carbon nanotubes where theory predicted a substantial spin-orbit 
coupling (a few hundred /leV) due to the curvature of the tube [78, 79, 80] which has 
been nicely confirmed in recent transport experiments on carbon nanotube quantum 
dots [81]. Since the surface of graphene is less curved than that of carbon nanotubes, 
the spin-orbit coupling in graphene - due to ripples - should still be rather weak (roughly 
ten times less than the spin-orbit coupling due to curvature in carbon nanotubes [79]). 
(ii) Carbon has two stable isotopes: ^^C and ^^C. The natural abundance is 99% ^^C 
and 1% ^'^C. Since ^^C has nuclear-spin and ^^C has nuclear-spin |, the electron spin 
of the qubit can only interact with 1% of the nuclei via hyperfine interaction. This ratio 
can even be further decreased because it is possible to artificially make ^^C-enriched 
graphene. 

4-2.1. Spin relaxation due to spin-orbit interaction and electron-phonon coupling in 
graphene quantum dots. In the long- wavelength limit, the spin-orbit interaction in 
graphene can be written in the following form [82, 83] 

Hsoi = ^iTCr^Ss + Aji{T(7^Sy - CTyS^), (17) 

where Sj denotes the Pauli matrix acting on the real spin, r = ± the valley degree 
of freedom, and is the Pauli matrix associated with the pseudospin due to the two 
sublattices of graphene. The first term (proportional to A/) has been coined intrinsic 
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spin-orbit interaction and is due to graphene's honeycomb lattice. The second term 
(proportional to A^) is of Rashba-type meaning that it is only present if the mirror 
symmetry about the plane is broken. This can be done by a perpendicular electric field, 
the interaction with the substrate, or ripples. 

The electron-phonon coupling in graphene is characterized by two mechanisms 
called deformation potential mechanism (DPM) - characterized by Vi below - and 
bond-length change mechanism (BLCM) - characterized by V2 below [84, 85, 86]. The 
DPM is diagonal in pseudo-spin space whereas the BLCM is off-diagonal in pseudo-spin 
space. Thus, the electron-phonon coupling Hamiltonian can be written as 



For spin relaxation in the low energy regime only acoustic phonons matter. There are 
three different acoustic phonons: (i) longitudinal acoustic (LA), (ii) transversal acoustic 
(TA), and (iii) transversal out-of-plane (ZA) phonons. Note that both LA and TA 
phonons have a linear dispersion whereas ZA phonons obey a quadratic dispersion. The 
specific form of Vi and V2 in Eq. (18) depends significantly on the type of phonons 
considered, see [80, 84, 85, 86, 87, 88, 89] for concrete examples both for graphene and 
carbon nanotubes. 

In a recent article, Struck and Burkard have shown that to lowest order in the spin 
orbit interaction as well as the electron-phonon coupling in flat graphene only the LA 
phonons matter for spin relaxation via the DPM while LA and TA phonons matter for 
spin relaxation via the BLCM [89]. They find spin relaxation times Ti between 10~^s 
and 10~^s depending on the apphed external magnetic field B that sets the level spacing 
between the two Kramers partners that form the qubit. This result should be compared 
to the corresponding analysis in GaAs [90, 91] as well as carbon nanotube quantum dots 
[80]. Interestingly, the predicted numbers for Ti are all in a similar range for different 
types of quantum dots but the B-dependence of the spin relaxation rate is totally 
different in different host materials. This is due to different symmetries and processes 
that dominate the relaxation. It confirms that a deep (microscopic) understanding of 
spin relaxation and dephasing in graphene quantum dots is crucial for an optimized 
qubit operation. Up to date, the spin decoherence time T2 due to spin-orbit interaction 
in combination with electron-phonon coupling has not been calculated for graphene 
quantum dots. For GaAs [91] and carbon nanotube quantum dots [80], it has been 
shown that T2 — 2Ti to leading order perturbation theory in spin-orbit interaction and 
electron-phonon coupling. 

All results of this section correspond to single layer graphene. In practice, it 
might turn out that spin qubits in bilaycr graphene quantum dots arc more stable than 
their single-layer counter parts. Rather little is known about the spin-orbit interaction 
[92, 93, 94] and electron-phonon coupling in bilayer graphene. Up to now, nothing is 
known about spin relaxation and dephasing in bilayer graphene spin qubits. We expect 
this to be a very active area of future research. 




(18) 
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4.2.2. Spin decoherence due to hyperfine interaction in graphene quantum dots. First 
investigations of hyperfine interaction in carbon nanosystems have been done on 
fullerenes - based on NMR theories and experiments [95] . Later on, ab initio calculations 
of the hyperfine interaction in small graphene nanoflakes have been carried out by 
Yazyev [96]. However, only very recently analytical calculations of the nuclear-spin 
interactions of electrons confined to graphene and carbon nanotube quantum dots 
have been presented [97]. The graphene part of the latter results will be reviewed 
in this section. Interestingly, pioneering experiments on nuclear-spin interaction and 
electron-spin dynamics in carbon nanotube quantum dots report an unexpectedly strong 
hyperfine interaction of order 100 /xeV in ^^C-enriched nanotubes [98, 99]. Unfortunately, 
similar data in graphene quantum dots does not exist up to now. Therefore, further 
experimental as well as theoretical work is needed to fully understand the role of nuclear- 
spin interaction on spin relaxation in carbon-based quantum dots. 

In principle, there are three terms that couple the spin of the confined electron 
to the nuclear spins: (i) the Fermi contact interaction, (ii) the anisotropic hyperfine 
interaction, and (iii) the coupling of electron orbital angular momentum to the nuclear 
spins. These interactions are represented by the Hamiltonians [100] 

hi = ^Y ^'^^^ ^^"'^ ^ ■ ^^^^ 
, fe /^o ^ ^ 3(nfc-S)(nfc-Ifc)-S-Ifc 

^'^^ r3(l + d/r,) ' (2°) 

respectively, where 7^ = 2/1^, 7^^. = gji^^N, I^b is the Bohr magneton, gj^ is the 
nuclear g- factor of isotopic species jk, fJ^N is the nuclear magneton, /xq is the vacuum 
permeability, = r — is the electron-spin position operator relative to the nucleus, 
d ~ Z X 1.5 X 10"^^ m is a length of nuclear dimensions, Z is the charge of the nucleus, 
and Uk — Tk/rk- S and = x p denote the spin and orbital angular-momentum 
operators (with respect to the k^^ nucleus) of the electron, respectively. The cutoff 
1 + d/rk avoids unphysical divergences from expectation values of the Hamiltonians /ig 
and /13. 

In order to determine an effective spin Hamiltonian, one first needs to calculate 
matrix elements of h^, h2, and h^ with respect to the electron states |\E'cr) = l^a] Ucr)\o-) 
where {r\^a]Ua) = $CT(r)ucr(r). Here, ^^(r) is the envelope function and Ua{r) the 
Bloch amplitude which has the periodicity and, thus, knows about the symmetry of 
the honeycomb lattice. In the following, the envelope function is not further specified 
and the Bloch amplitude is approximated by hydrogenic orbitals, see Ref. [97] for more 
details. In flat graphene, the electron states are p-type, whereas in nanotubes and curved 
graphene they are sp-hybridized. We restrict ourselves to flat graphene in the remainder 
of this section. Then, the matrix elements with respect to h^ and h^ vanish, and only 
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cr\h2\^ a') is finite. It gives rise to the following effective spin Hamiltonian 

HuFi = Yl (^^^'4' + ^l^'^k + KS'Il) (22) 

k 

with coupling constants A^^ — Aji>o|$o-(rk)p (where Vq is the area of a primitive unit 
cell). Interestingly, the hyperfine interaction is asymmetric in space with [97] 

where the z direction is perpendicular to the graphene plane. In Eq. (23), Z^s labels 
the effective screened nuclear charge and ao is the Bohr radius. It has recently been 
shown, in the context of transport calculations on carbon-based quantum dots in the 
so-called Pauli blockade regime, that the hyperfine Hamiltonian gives rise to a coupling 
between the valley degree of freedom of the electron and the nuclear spins [101]. Such a 
coupling is not important for the interaction between the electron spin and the nuclear 
spins discussed here. 

We now address the spin decoherence properties of a single electron spin in a 
graphene quantum dot coupled to a bath of nuclear spins in the situation where a 
"large" magnetic field B.. is applied perpendicular to the graphene plane. Specifically, 
this means that b = g^sBz ^ Aj which in fact corresponds to a rather moderate 
magnetic field of about 5 mT. Under these circumstances, the main source of decoherence 
is pure dephasing due to nuclear-field fiuctuations along the spin quantization axis in 
2;-direction governed by the Hamiltonian 

H = (h + K)S' (24) 

with hz = ^fc-^fc- The total number of nuclei can be written as = + A^is where 
A^i2 is the number of ^^C nuclei and A^i3 is the number of ^^C nuclei. For Ni^ ^ 1 we 
can use the central limit theorem (compare with Ref. [102]) which yields the following 
Gaussian dynamics for the transverse electron spin 

(5+), = (5+)oe-*'/^c, (25) 

where 5'"'" = ^S"^ -|- iS^. The latter equation is written in the rotating frame with 
rotation frequency {u + pr]Az/2)/h, where u — h — with the nuclear Zeeman energy 
— QnI^'mBz- The characteristic electron-spin decoherence time is given by [97] 
~ - 2/1 N 

where the polarization p of the nuclear-spin system (0 < p < 1) is determined only by 
the distribution of spin-up and spin-down ^^C nuclei. In Fig. 10, we plot as a function 
of the relative ^^C abundance showing very promising decoherence times for graphene 
quantum dots. The number for the decoherence time for the natural ^^C abundance 
fc > 80//.S [97] should be compared to the predicted 5ns for GaAs [102] in the absence 
of nuclear-spin polarization. This is a gain of more that 4 orders of magnitude. 
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Figure 10. Electron-spin decoherence time fc in graphene. fc is plotted as 
a function of the relative ^■^C abundance Ni^/N, under the condition of a magnetic 
field Bz ^ 5 niT perpendicular to the graphene plane. We assume no nuclear-spin 
polarization and estimate the total number of nuclei to be N — 4 x 10^. Inset: 
Electron-spin dynamics in graphene with ^"^C abundances of 1% (solid blue curve) 
and 99% (dashed red curve), in the rotating frame (after Ref. [97]). 



5. Conclusions and Outlook 

We have reviewed the efforts towards the creation of quantum dots in graphene with 
the emphasis of using them as a host for spin qubits. Graphene has great potential 
for being an ideal candidate for spin qubits due to its low intrinsic spin-orbit coupling 
and the sparse amount of nuclear spins. We discussed theoretically in detail bound 
states in gate-tunable graphene quantum dots realized in graphene nanoribbons and 
in gapped single-layer and bilayer graphene. In contrast to quantum dots realized in 
edged graphene flakes, gate-tunable quantum dots are defined electrostatically and not 
by the physical edge of a graphene sample. This allows to controllably break the valley 
degeneracy, a prerequisite for spin-based quantum computing, e.g. by using a magnetic 
field. We have also discussed quantum manipulation of spin qubits in such dots as 
well as very recent theoretical studies on the consequences of spin-orbit interaction 
and hyperfine interaction with nuclei for spin-relaxation and spin-decoherence. Since 
theoretical as well as experimental efforts have so far focused mainly on single-layer 
graphene quantum dots, we believe that a great deal of interesting work will be devoted 
to bilayer graphene. Bilayer graphene is potentially superior to single-layer graphene due 
to the creation of a tunable bandgap by electric fields which allows for an all electrical 
control of graphene quantum dots. 
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